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Abstract 



We derive the spectral representations of QED 3-point functions and then 
explicitly calculate the 3-point spectral densities in hard thermal loop approx- 
imation within the real time formalism. The Ward identities obeyed by the 
retarded and advanced 2- and 3-point functions are discussed. We compare 
our results with those for hot QCD . 
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I. INTRODUCTION 



Spectral densities are important quantities in finite temperature field theory |],|2[]. But 
they are not easy quantities to evaluate perturbatively at nonzero temperature, especially for 
many-point spectral densities. The two-point spectral densities have been widely studied and 
applied to QGP studies In Ref. J/J the Cutkosky rules for calculating the imaginary 

parts of thermal Green functions using the formalism of Thermo-Field Dynamics (TFD) 
was presented. Recently these cutting rules were reexamined within the Closed Time Path 
(CTP) formalism [rU,|TT] and given a simple physical interpretation [O. The imaginary 
time cutting rules for calculating two-point spectral densities were investigated in || . Three- 
point spectral densities, on the other hand, have so far received less quantitative attention. 
The only available calculation for gauge theories is published in |13| where the 3-point 
spectral densities for pure gluon dynamics were calculated in the Hard Thermal Loop (HTL) 
approximation using the imaginary time formalism (ITF) |L4j]. The spectral representation 



of three-point functions for selfinteracting scalar fields were discussed in Refs. 1 15-18 



Hard Thermal Loops (HTLs) are gauge invariant and satisfy simple abelian Ward identi 
ties [[14]]. These remarkable properties have triggered many interesting investigations |T9|-|24 



The computation oh HTLs is generally fairly technical because of their complicated momen 
turn and energy dependence, but it can be simplified by using the Ward identities HI^PS" 



All HTLs can be derived from a generating functional based on an effective Lagrangian 
19| . They describe classical aspects of hot field theories and can thus also be obtained from 



classical kinetic equations |26|j23|j24]| . However, as Taylor stated in [13], for many purposes 
one does not need the HTL amplitudes themselves, but only their discontinuities which are 
described by spectral densities. These spectral densities also provide a natural connection 
between the Green functions in the real time formulation (RTF) of thermal field theory and 
the ITF. While the HTL resummation method was developed within the imaginary time 
formalism (ITF), with primary attention focussed on equilibrium properties of hot field 
theories, realistic physical systems are frequently out of thermal equilibrium and require 
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calculations in real time. Recently this has motivated increased interest in the real time 
formulation of thermal field theories and its non-equilibrium extensions. 

In the present paper we therefore study the spectral functions and Ward identities for 
finite temperature QED in the real time formalism (RTF). Extending the recently derived, 
completely general spectral representation of the real-time 3-point vertex function at finite 
temperature from scalar field theory to the case of QED, we then show how to implement the 
HTL approximation in real time. As shown in [p7H , for 2-point functions this procedure al- 
lows for a simple generalization of the HTL resummation scheme to general non-equilibrium 
situations. Within the HTL approximation we evaluate explicitly the 3-point spectral den- 
sities and derive a set of finite temperature Ward identities between the real-time 2- and 
3-point functions. Although one needs in general two independent spectral densities to de- 
scribe the real-time 3-point vertex at finite temperature |I7|,|18|, the two are shown to become 



degenerate in the HTL approximation. This agrees with previous findings in Refs. (|T3| , |23 
Our Ward identities between the real-time Hard Thermal Loops in QED also agree with the 
general real-time finite temperature Ward identities recently derived in Ref. 



Throughout this paper we will use the CTP formalism [ItJ in the form given in 



Refs. JTT|,pT|. In this representation of the real-time formalism the bosonic single-particle 



propagator in momentum space has the form 



D n D 12 

Dip) = ( ) (1) 

D 2 i D 22 



with 



iD n (p) = (iD 22 )* = iV bf^J + ("(Po) + I) p{p) , (2a) 

i D 12 (p) = n(po) pip) , (2b) 

iD 21 ip) = (l + nip ))pip). (2c) 

Here n(po) is the thermal Bose-Einstein distribution, 



and p(p) is the two-point spectral density which for free particles is given by 

pip) = 2tt sgn(p ) 5{p 2 -m 2 ). (4) 
The fermionic 2x2 propagators are given by 

Sab{p) = {p' + m)S ab , (a, b = 1,2), (5) 



with 



• &i(p) = (i S 22 y = iV (^4^) + + I) PiP) > 

i Sn(p) = -"(Po) p(p) , (6b) 

i5 2 i(p) = (l-n(po))p(p). (6c) 



Here n(po) is the thermal Fermi-Dirac distribution 

1 

e Ppo + i 



"(po) = , 1 ' ( 7 ) 



and for free fermions p(p) is again given by (H). 

The paper is organized as follows. In Sec. [TJ we derive the spectral representations for 
the retarded 3-point functions in QED. In Sec. |TT1| we evaluate the 3-point spectral densities 
for QED in the HTL approximation. In Sec. |V| we derive the RTF Ward identities between 
the 2- and 3-point HTL amplitudes in QED. A short summary is given in Sec. 0. 

II. SPECTRAL REPRESENTATION OF THE 3-POINT VERTEX IN QED 

In this Section we shortly review some useful relations among the different thermal 
components of the QED 3-point functions and derive their spectral representation. Similar 
relations for the 3-point vertex in 3 theory have been reported in the literature [[7|,|l5|,[n| in 
different notation. Our procedure here follows the notation developed in |18j for 3 theory. 
We consider the 3-point vertex function in QED shown in Fig. 1. The three incoming 
external momenta are k\ = p, k 2 = q, and k 3 = —p — q. 



Fig. 1: 3-point vertex in QED 



Within the real time formalism, the truncated three-point function for QED, G^ bc (x,y, z) 
has eight thermal components which satisfy |7|,|T0"1 



Gabc 

a,b,c=l 



0. 



Using the KMS condition one finds in momentum space [ [30] 



^2) ^3) — —Gili{k\ } k 2 , ks) — G 222 (ki, k 2 , k 3 



G^ 2 i{ki,k 2 ,k 3 ) 
G% 11 (k 1 ,k 2 ,k 3 ) 



-G^ih, k 2 , k 3 ) = -e^ G^ih, k 2 , k 3 ) . 
-£2*1 (&i> ^2, h) = -e^ 1 Gl 22 (kt, k 2 , k 3 ) . 
-G M *112(/ci, fc 2s fc 3 ) = e^ 3 (^221(^1 j fc 2 , fcs) 



(9a) 
(9b) 
(9c) 
(9d) 



where G M represents "tilde conjugation" of G M (see PI). Note that in Eqs. |)b|,|9c]) the 
last equation involves an additional minus sign relative to the scalar case ||18|| , due to the 
fermionic legs. 

One can construct "retarded" vertex functions from the above eight thermal components 
according to 



rip r<P 1 rip 1 rP 1 r<P 

— ^111 t LT112 ~r LT211 "r "^212) 



r<P rP 1 rP 1 rip _i_ ^ 

^ffi — ^111 + ( - T 112 + ^121 + 



P 

122) 



LT i?,o — Wll 



rip 1 /^M 1 /^A* 
^121 t LT211 t <J"221) 



(10a) 
(10b) 
(10c) 



where is the vertex function which in coordinate space is retarded with respect to xq, 
is retarded with respect to z , and is retarded with respect to y . The legs linked to 
x, y, z correspond to an electron, electron, and photon, respectively, and the corresponding 
distribution functions are hi = h(ui), h 2 = n(u; 2 ), and n 3 = n(u 3 ). Inversion of Eqs. ( |I0D 
with the help of (|8|) (^) yields expressions for the thermal components G^ bc in terms the 
above retarded functions: 

a^MM = ^ y (ij (-J - + . s) (_ 1 i ) (, _\) (_\) 

+ <* (_\) G _\) (r;:J - \™ + *> G y W (_\) 
+H\)G- 2 J(- 1 1 ) - H^^W (-\)G?J. <»> 

with Aq = 1 — 2hi, N2 = 1 — 2fi 2 , and A3 = l + 2n 3 . The structure of this equation is similar 



to Eq. (35) in Ref. |fT7 ] ; the sign differences arise from additional minus signs in front of the 
Fermi distributions and from the fact that Eq. (pTTJ ) refers to the truncated vertex rather 
than the connected vertex studied in |Ty] . 



Following the same procedure as in Appendix A2 of Ref. [181] one derives the following 



spectral integral representation for the retarded 3-point functions in QED: 

G^(ui,u 2 ,u 3 ) = — / — ■ — - + - , (12a) 

r»< , _ -i r d ^ d ^ ( P2 P1-P2 \ , 19M 

G m {u Xl u> 2 , u 3 ) - — 2 J ^ Ui _ Qi+ it y~ _ Q2 _ it - ^ Z Qs _ ie ) , (12b) 

G% ) (u 1 ,u 2 ,u 3 ) = / - - + p: . 12c 

z7r J-00 W3 — il 3 + ie \u 1 — ill — ie oj 2 — "2 — J 

Here oj\ + 0J2 + uj 3 = fi x + Q 2 + ^3 = and the spatial momenta pi, p 2 , p 3 = — (pi— P2) 
are the same on both sides and have therefore been suppressed. The spectral densities are 
given by the following thermal components of the 3-point vertex in momentum space: 

p? = Im(G? 22 -^ n ), (13a) 
p^ = Im(G^ 12 -G^ 21 ). (13b) 
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One notes that the spectral representations for the truncated QED three-point functions 
have the same form as in scalar 3 theory |)TB| , except for the additional vector index. 



III. EVALUATION OF THE SPECTRAL DENSITIES IN HTL APPROXIMATION 

In this Section we calculate the three-point spectral functions pi , p 2 for QED in the hard 
thermal loop approximation | 14j| . Since HTLs are gauge invariant we can choose Feynman 



gauge for simplicity. 

From Eqs. (§) and ([13]) we have 

pi = Im (G? 22 - eft*G&) = -^rlm G? 22 . (14) 

n(po) 

For pi we thus must evaluate only the single Feynman diagram in Fig. 1 for a = 1, b — c = 2. 
Using standard real-time Feynman rules |T(J and the photon propagator in Feynman gauge 



one gets 



G? 22 (p,q,-p-q) = {-ig){igY ) j^iiDMs^jaiiS^s + q^liSnis-p)}^. (15) 
Inserting the thermal free propagators ©-(0) and extracting the imaginary part one finds 



ImG? 22 = -0*1—1- [ la ^-0 + m )YU+ 4+m)l a ] 



(2tt)< 

x 5(s 2 ) S((s + q) 2 - m 2 ) 8((s - p) 2 - m 2 ) sgn(s ) sgn(s + q ) sgn(s - p ) 
x n(s ) (| - n( s o + 9o)) (1 - H s o ~ Po)) (16) 

If the coupling constant g is small and the external momenta are soft, p,q ~ gT, and 
the electron bare mass m is much smaller than the temperature, the leading contributions 
comes from the hard loop momenta s ~ T |4]Jlil . For these we can thus neglect the external 



momenta and the mass m in the terms between square brackets, approximating them by 
7o^7/i^7a = — 4s M / + 27 M s 2 . After performing the integration over s° with the help of the 
function 5(s 2 ) = [5(s — s) + 8(s + s)]/2s, where s = Vs^, one finds 
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flip, Q, ~V-q) = -rr(^(P, l) + S m (p, ?)) - (17a) 
n(p ) v ' 

/c?^s 1 
(27r) 2i Sgn ^ + q °) Sgn ^ ~~ Po) 4s ^AsO=s 

x5((s + g ) 2 - E 2 s+q ) 5{{s - p ) 2 - E 2 _ p ) 

xn(E s ) (i - fi(s + g )) (1 - n{s - Po)) , (17b) 
f d n ~ l s 1 

S^(p, q) = - J ^ 27r )n-3 2i sgn ^ s ~ sgn ( s + Po)4s M /| s o = _ 5 
x5((s - g ) 2 - E 2 +q ) 5((s + po) 2 - El p ) 
x (n(-s)) (| - n(-s + g )) n(s + p ) . (17c) 



Here E s+q = Jm 2 + (s + g) 2 , -E s _ p = ym 2 + (s — p) 2 , and in ( |17c| ) we used the identity 
n(— x) + n(x) = —77, 77 = ±1 for bosons and fermions. 

To simplify the notation it is convenient to introduce the 4-vectors V± = (l,±|). For 
hard loop momenta the arguments of the 5-functions in ( |T7|) can then be written as 

(s + q) 2 n ^ ±2s q ■ V± = ±2s (q°T\q\ cos 6') , (18a) 



ps ±2sq 


v± 


= ±2s 




\q\ 


cos 9') 




v± 


= T2s 






cos 



(s-pY , PS ^2sp- 14 = T2s p u T Iplcosfl , (18b) 

and we can replace sgn(s ±p°) « 1 ^ sgn(s ± g°) under the integral. 9 and 0' are the angles 
between s and p, q, respectively. These approximations decouple the angular and radial 



integrations WM; we obtain 



A"{p,q)^a{p )u IM {p,q), (19a) 

B»(p,q)&b(p )co»(p,q), (19b) 

1 f°° 4s 4 / \ 

a ^ = 2~nlo rf ^(2i)3 n ^H^ -«(* + ?o)J(l -n(s-po)), (19c) 

1 f 00 4s 4 / \ 

& (Po) = -7^ y o tis^z^n(-s) [2 -n(-s + g )Jn(s + p )), (19d) 

w"(p,g) = JdnV^f + 5(q- V+)5(p- V+) = J dttV^f _ 5(q • VL) 5(p ■ V-) . (19e) 



Except for the additional spinor structure the angular integral (|19ej ) is identical with the 



one found by Taylor [13] for the spectral density of the 3-gluon vertex in hot QCD. 
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The integrands in Eqs. ( |19c| , |19d|) contain up to three powers of thermal distribution 
functions. By using following identities 



n(uji)n{ijj 2 ) = n(a>i + u; 2 )(l + n{uj x ) - h(u 2 )) , 
n(uji)h{ijj 2 ) = n(u)\ + w 2 )(l - n(o»i) - n(u> 2 )), 

one shows that the cubic terms disappear and that a and 6 reduce to 

o(po) = ^ y o dss(|n(p ) [n(s) + n(s - po)] 

-n(po + gb)[n(-go)[n(s) +™0 + g )] -n(po)[n(s) +n(s-po)] 

1 /"°° / 
KPo) = ^ y o dss(yi,(p ) [n(s) + n(s+p )] 

-n(p + <?o) n(-?o) + n(s - go)] - n(po) K«) + n(s + p )] 



(20a) 
(20b) 



(21a) 



(21b) 



The remaining integrands are linear in the s-dependent thermal distribution functions 



31,18 



After evaluating the integral over s in the limit Po/T <C 1, q /T 1 we obtain to leading 
order in the coupling constant g 

,3 



pi(p,e) 



g 



n(po) 
gm 2 p7r 



(A»(p,q)+B»(p,q)) 

J dnV^f + 5{q-V + )5{p-V + ) 
J dVt V^Y _ 5{q ■ V.) 5{p ■ V-) 



[22) 



Here mp = ^ is the thermal electron mass. 

The other spectral density p 2 is determined from 

1 



P2 



Im G£ 12 . 



(23) 



By inspection of the corresponding labelling of the diagram in Fig. 1 one observes that 
G 2 i 2 (p, q, —p — q) is obtained from Gi 22 (p, q, —p — q) by exchanging the two electron legs 
with the external momenta p and q and routing the internal momentum s in the opposite 
direction. The resulting loop integral then becomes identical to the one before, Eq. (|T5|), 
and we find 



P2 (P,Q) « Pi (<?, P) ~ Ph TL (p, ?) , (24a) 

2 

p£ tl (p, g) « ^ / dn v£ f + 5(p-F + ) %-v + ) . (24b) 



This agrees with the observation by Taylor within [T3[ that in QCD in HTL approximation 



the two independent spectral densities for the 3-gluon vertex degenerate. 
Considering the ^-functions in Eq. ( P4[) one easily shows that 

P^Ph TL (P, q) = = ^Ph TL (p, ?) (25) 
i.e. in HTL approximation the QED spectral density is transverse with respect to the 



external momenta. This nice feature was also noted in [13] for hot QCD. 



IV. REAL-TIME WARD INDENTITIES AMONG HARD THERMAL LOOPS 

Let us insert the above spectral densities into the spectral representations for the retarded 
3-point functions: 

m 2 r V 11 V 

47r J (V + ■ ki — ie){V + ■ k 2 + ie) 

G&iuu wa, a* = ~W^ dQ + ; + , (26b) 

47r J (V+ ■ fci + «e)(V+ ■ k 2 — ie) 

(V+ ■ h - ie)(V+ ■ k 2 - ie) 

The QCD-analogue of last two expressions were previously obtained by Blaizot and Jancu 
23fl from a set of classical kinetic equations. If the external momenta are of order gT, 



G£> 1; u 2 , u 3 ) =ig~^ dn—- . (26c) 

4:71 J (V-L. ■ Ki — ie (V4- ■ ko — ie) 



Pi ~ P2 ~ power counting reveals that all three retarded 3-point functions are of order 
g, i.e. of the same order as tree vertex. The HTLs fl26"D for the QED 3-point vertex in the real 
time formalism thus require resummation within perturbation theory, like the corresponding 
HTLs in the imaginary time formalism |14,|5[. 



In the next step we derive the real-time analogue of the well-known Ward identities 



between the HTL amplitudes in the ITF formalism [25|,19|. Due to the matrix structure of 



the real-time thermal Green functions, the Ward identities also become matrix equations. 
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The zero temperature Ward identities for QED can be written both in differential and 



in integrated form [29 



^ = - g S{v) r M (p, -P, 0) S(p) , (27) 



or 



9 W) - S(p)) = i S(p) i(p' - p)%(p, - V \V - p)\ S(p) ■ (28) 

Multiplying by S~ 1 (p) from the left and by S^ 1 (p') from the right and using the Schwinger- 
Dyson equation 

S-^p) = ^-m-E(p), (29) 
the Ward identity takes the form 

(p - pYG^p, -p',p' -p) = -ig [Eft/) - E(p)] , (30) 
where is the vertex correction from loop diagrams, defined by G M = V ^ — T Q , with 

At finite temperature the propagators and vertices have matrix structure. The general- 



ization of (30) to finite temperature in the real time formalism was recently performed in 
Ref. [pSfl . The authors of this paper used a different version of the real time formalism and 
considered the vertex with the photon attached to the ingoing leg; translated into the CTP 
framework and for the vertex with the photon on the outgoing leg their result reads: 

(p' - p)tiG^ bc (p, -p',p' -p) = -ig [5 bc Y>ac(p') - 6 ac J: cb (p)} , (31) 

where a,b,c = 1, 2. This can be rewritten in terms of the retarded and advanced amplitudes; 
the result agrees with what we find for the HTL amplitudes below. 

In order to examine the relations among the 2- and 3-point HTLs in the real time 
formalism, we calculate the retarded electron self-energy JlTJ: 



. . 9 r d 3 k r , , , . . . , . . . / n(u') + n(u) \ , , 

MP )=9 2 duduj'p f (uj,p + k)p B (u',k ' 32 

J (2,71) J \po + uj ' — LU + iel 
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where pp and pB are the two-point spectral densities of the fermions and bosons, respectively. 
Inserting the free particle spectral densities 

p B (u', k) = sgn(a/) 5(uj' 2 - El) , 
p F (u,p + k) = (^7° + k' • 7 + m) sgn(cu) <^ 2 - E 2 p+k ) , (33) 

with Ek = Vk 2 , E' k = \fk' 2 + m 2 , k' = k + p, and evaluating the integral at in the HTL 



approximation and using dimensional regularization p2| , we obtain 

E R (p) = a(p)^ + b(p) lo , (34) 

where 

a(p ,p) = ^(l-^Jdn 



V + ■ p + ie 



b(po,p) = ^f(po + f[dn 1 ) . 

p z \ AlT\ J V + ■ p + iej 



(35) 



Here V + = (1, V),V = j^, and the integration is over the direction of the unit vector V. 
Some further algebra then leads directly to 

z R (p) = Hi f d n y+ , . (36) 

y > 4tt J V + -p + i6 y ' 

The advanced electron self-energy in HTL approximation is computed similarly as 

Y 



Sa(p) = ^ fdSl V . ■ (37) 
An J V + ■ p — ie 



From this one obtains 



47r 7 (V+ ■ p + «e)(V+ • p' + ze) 

47T j (V+ ■ p — «e)(V+ ■ p + ie) 

S A (p) -E A (p) = --^ / dO— —77 y—r- 38c 

47r J (V+ ■ p — ie){V+ ■ p — ie) 



If we set in ( p6|) fci = p, &2 = — p', k% = p' — p and then compare it with Eqs. (|38| ) we find 
the following relations between the hard thermal loop contributions to the 2- and 3-point 
functions: 
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(p - p)nGm(p, -P',P' -P) = -w\Zr{p') ~ ^r(p)} , (39a) 
(P ~ p)»Gr(p, -p',p' -P) = -w[Za(p') ~ ^a(p)} , (39b) 
tf-pTG&b, -p',p'-p) = -ig\EA&) - E*(p)] . (39c) 



These relations are structurally similar to the zero temperature Ward identity (|5D|) and agree 
(up to notational differences) with the general finite temperature Ward identities in RTF 
given in Eq. (3.17) of Ref. ||. 



V. CONCLUSIONS 

We studied the QED 3-point vertex function at finite temperature in the CTP real-time 
formalism. This formalism has recently gained increased popularity because it allows for 
a generalization to non-equilibrium situations as encountered, e.g., in the initial stages of 
heavy-ion collisions, and it avoids the need for analytical continuation which plagues the 
imaginary time formalism. We started by giving a set of useful relations among the eight 
thermal components of the real-time vertex function. We then derived spectral integral 
representations for the three retarded 3-point functions and calculated the corresponding 
spectral densities explicitly at 1-loop order in the HTL approximation. In this approximation 
the two independent spectral densities become degenerate and are transverse with respect 
to all three external momenta. 

Inserting these HTL spectral densities into the spectral representation we obtained three 
retarded 3-point vertex functions, two of which turned out to be identical to the corre- 



sponding QCD-analogues derived in [23| from a set of classical kinetic equations in the 
long- wavelength limit. By contracting the three 3-point vertex functions with the momen- 
tum vector of the photon, we obtained a result which we could compare with the fermion 
HTL self energy. The result was a set of real-time Ward-Takahashi identities at finite tem- 
perature which generalize the zero temperature Ward identity and agree with the recently 



derived finite temperature identities of Ref. [EB . 



13 



Due to the matrix structure of the real-time thermal Green functions, there is a whole 
class of finite temperature Ward identities which relate retarded and advanced vertex func- 
tions to combinations of retarded and advanced fermion self energies. As first observed by 
D'Olivo et al. pSfl , if the ingoing (outgoing) fermion leg has the largest time, the Ward 
identity involves only the retarded (advanced) fermion self energies; if the photon leg has 
the largest time, both retarded and advanced fermion self energies are involved. 

Clearly, no fundamentally new physical results were derived in this paper; this was not 
our goal. What we have achieved is a consistent real-time representation of the 2- and 3-point 
functions in finite temperature field theory, taking into account the full matrix structure aris- 
ing from the doubling of degrees of freedom, which makes their analytic structure explicit 
and can thus serve as a basis for non-perturbative resummation schemes at finite tempera- 
ture. We have tested the consistency of the formalism in the context of HTL resummation 
within thermal equilibrium QED. We expect the formalism to provide a stable basis for an 
extension to the real-time dynamics of non-equilibrium systems. 
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